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ABSTRACT 



In recent nume r ical simulations (jMatsubayashi et al.l 120071 ; 



Lockmann. fc Baumgardtl 120081 ). it has been found that the eccentricity of 
supermassive black hole(SMBH) - intermediate black hole(IMBH) binaries 
grows toward unity through interactions with stellar background. This in- 
crease of eccentricity reduces the merging timescale of the binary through the 
gravitational radiation to the value well below the Hubble Time. It also gives 
the theoret ical explanation of the e xistence of eccentric binary such as that 
in OJ287 dLehto. k Valtonenl Il996l : IValtonen et~aD l2008h . In self-consistent 
N-body simulations, this increase of eccentricity is always observed. On the 
other han d, the result o f scattering experiment between SMBH binaries and 
field stars ( IQuinlanlll996l ) indicated no increase of eccentricity. This discrepancy 
leaves the high eccentricity of the SMBH binaries in iV-body simulations 
unexplained. Here we present a stellar-dynamical mechanism that drives the 
increase of the eccentricity of an SMBH binary with large mass ratio. There 
are two key processes involved. The first one is the Kozai mechanism under 
non-axisymmetric potential, which effectively randomizes the angular momenta 
of surrounding stars. The other is the selective ejection of stars with prograde 
orbits. Through these two mechanisms, field stars extract the orbital angular 
momentum of the SMBH binary. Our proposed mechanism causes the increase 
in the eccentricity of most of SMBH binaries, resulting in the rapid merger 
through gravitational wave radiation. Our result has given a definite solution to 
the "last-parsec problem" 
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1. Introduction 

Most of galaxies contain SMBHs at their 
centers. If a galaxy collides with another 
galaxy, the merging remnant hosts two 
SMBHs and th e y form a binary system. 
Begelman et al.l (ll980l ) discussed the pos- 
sibility of the formation of SMBH binaries 
and their evolutions. They concluded that 
the evolution of the semi-major axis stalls 
after they ejected out all the nearby stars 
(loss cone depletion) and the binary is un- 
likely to merge within the Hubble time, 
as far as we consider the stellar dynamical 
effects and the gravitational wave radia- 
tion only. This problem is called as "last 
parsec problem" since the typical separa- 
tion of a stalled binary is around one par- 
sec. Many researchers studied this prob- 
lem and confirme d that the separation of 
the binary stall (IMakino fc Funatol 12004] ; 



Berczik et al.l 120051 ). Most of previous 



works on the evolution of the SMBH bi- 
nary focused on the e volution of the semi- 



major axis. Recently, iMatsubayashi et al. 



(120071 ) performed iV-body simulations of 
SMBH-IMBH binaries with the mass ra- 
tio of 1000:1, and found that the eccen- 
tricity of the binary grows toward unity 
after the loss cone depletion. Because 
of this increase of eccentricity, the merg- 
ing timescale through the gravitational 
radiation goes below 10 7 years even if 
the loss cone is depleted. Similar result 
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was ob tained by lLockmann. fc Baumgardt 
( 120081 ). and it has been reported that un- 
der certain circumstance SMBH binaries 
of mass ratios close to 1:1 show similar in- 



crease of the eccentr icity (lAarsethl 12008 



Berentzen et all 120091 ). 



However, this increase of the eccentric- 
ity is difficult to understand, because the 
measured differential cross sections of in- 
teraction between an SMBH binary and 
field stars indicated that eccentricity of a 
SMBH binary would not i ncrease throug h 
interactio n with field stars (|Quinlanlll996l) . 
Recently, ISesana et al.l ( 120081 ) pointed out 
that the bound stars around SMBH binary 
play important role in the evolution of the 
eccentricity but did not give the explana- 
tion about the mechanism of the eccentric- 
ity growth. 

If the increase of eccentricity actually 
occurs in real SMBH binaries, it nat- 
urally explains the o rbit of the binary 
in the quasar OJ287 dLehto. fc Valtonen 



1996; Valtonen et al. 20081). The central 



SMBH of the quasar OJ287 has the mass 
of 1.8 x 1O 1O M , and the secondary black 
hole has the mass ~ 1O 8 M , orbital period 
of 12 yrs, and eccentricity of 0.66. Cur- 
rently, this SMBH binary is in the process 
of circularization due to the gravitational 
wave radiation. Thus, one would expect 
that OJ287 had much larger semi-major 
axis and eccentricity very close to unity in 
the past. 

In this letter, we describe the mecha- 
nism which drives the increase of the ec- 
centricity. In section 2, we describe the nu- 
merical experiment we performed. In sec- 
tion 3, we present the interpretation of the 
numerical result. We give the discussion 
and summary in section 4. 
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2. Numerical Simulation 

In order to clarify the mechanism which 
drives the growth of the eccentricity of 
unequal-mass SMBH binaries, we per- 
formed iV-body simulations and investi- 
gated the change in orbital property of field 
stars due to interaction with the SMBH bi- 
nary. Our initial mo d el is s imilar to that of 
Matsubayashi et al.l ( 120071 ). Galaxy den- 
sity profile is given by 



p(r) 



n 



M n rl 



4tt r 3_r ?(r 2 + r 2 ) 7 ?/ 2 " 1 



(1) 



where r$ is the scale length and is 
the total mass of field particles, and we 
set r = 60pc and M v = 1.25 x 10 9 M Q . 
From observations, the slope of central den- 
sity of a gian t elliptical galaxy is shal- 
lower t han —1 (Lauer et al. 19951 ; Faber et al 



19971 : iKormendv et all 12009^ . This fea 



ture is well understo od by numerica l sim- 



ulations and theory (lEbisuzaki et al.l 11991 



Nakano fe: Makino I Il999al lbl). In our simula- 
tion, the density is p oc r -0,75 cusp. The ve- 
locity profile is isotropic and consistent with 
the density profile. The number of field stars 
./V are 16384 and 32768. We performed five 
runs for each N. The heavier SMBH parti- 
cle whose mass is M p = 1O 1O M is placed at 
the center of the galaxy, and the lighter one 
whose mass is M s = 1O 8 M is placed m an 
eccentric orbit with eccentricity e of 0.2 and 
semi-major axis a of 20pc. At this radius the 
enclosed stellar mass is comparable to that of 
the secondary mass and the stall of the evo- 
lution of semi- major axis takes place quickly. 
The mass of field particles 3.8 x 10 4 M Q in the 
case of N = 32768. Our simulation roughly 
mimics the OJ287 system. 

The time integration sche me we used is the 
fourth-order Hermite scheme (Makino &: Aarseth 



19921 ) for interactions between the field-field 
particles and the secondary SMBH - field par- 
ticles. The forces from the primary SMBH 



are inte grated with the sixth-ord er Hermite 
scheme ( Nitador Sz Makino 20081 ) to reduce 
the integration error. To save t he calculation 
time, we used GRAPE-DR (jMakino et al 



2007) for calculating interactions between the 
field- field particles. 

Figure 1 shows the evolution of the semi- 
major axis a, the eccentricity e and the merg- 
ing timescale due to the gravitational w ave ra- 
diatio n Tqw of a SMBH binary using Petersl 
formula, for an A = 32768 run. The 
evolution of the semi-major axis slows down 
and deviates from the theoretical prediction. 
At the same time, the eccentricity starts to 
grow, and reaches ~ 0.99 by T = 260 Myrs. 
This increase of eccentricity causes the dras- 
tic decrease of the merger timescale through 
gravitational wave radiation. In this case, at 
T = 180 Myrs, the merger timescale is 10 
Gyrs, but by T = 250 Myrs the timescale 
shrinks to 10 Myrs. We did not include the ef- 
fect of the gravitational wave radiation, but if 
we did, the SMBH binary would have merged 
very quickly. 

Figure 2 shows the evolutions of the ec- 
centricity for all runs. The evolutions in all 
runs inferred from Figure 2 are basically simi- 
lar. The variation of the results of N = 32768 
runs is smaller than that of N = 16384 runs. 
These appear to be no statistically significant 
difference between the results for N = 16384 
and 32768. 

3. Why does the eccentricity grow? 

We have carefully analyzed the behavior of 
the stars around the SMBH binary in our nu- 
merical simulation, and found that the follow- 
ing three-stage mechanism drives the increase 
in the eccentricity of the secondary. 

1. The presence of the secondary gives 
small, non-axisymmetric perturbation 
to the orbits of the field stars, which 
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would be purely Keplerian otherwise. 
Because of this perturbation the orbital 
motions of the stars become chaotic, 
and their angular momenta change 
chaotically. 

2. Since the orbital angular momenta of 
stars are not conserved, they can ex- 
perience close encounters with the sec- 
ondary, some of which result in the ejec- 
tion of stars. When a star is on a pro- 
grade orbit relative to the orbit of the 
secondary, it has a much larger chance 
of being ejected than when it is on a 
retrograde orbit. 

3. Since the secondary tends to scatter 
stars on prograde orbits, there will be 
more stars on retrograde orbits than 
on prograde orbits around the binary. 
However, the chaotic changes of angu- 
lar momenta of stars reduce the differ- 
ence between the number of stars on 
prograde and retrograde orbits. This 
means the angular momentum of the 
secondary is transferred to field stars 
and the eccentricity of the secondary in- 
creases. 

Note that in this process the transfer of 
angular momentum occurs between the bi- 
nary and field stars, while field stars are 
bound to the system, through secular pertur- 
bation. This secular perturbation is neglected 
in previous theoretical models on the behay - 
ior of SMBH binaries ( Begelman et al. 1980l ) 
or mode l s bas ed on scattering experiment 
(|Quinlan|[l~996h . and can give t he reasonable 



i nterp retation of the results of ISesana et al 
(|2f)08h . 



Figure 3 shows the time evolution of the 
specific angular momentum of a star that 
comes close to the secondary SMBH. The or- 
bital plane of the secondary SMBH is the 
x — y plane. One can see that none of 



the three components of the angular mo- 
mentum are conserved and the changes are 
complicated. These changes are due to 
the secular perturbation of the secondary 
SMBH. If the secondary SMBH had a circu- 
lar orbit, the secular perturbation would not 
have any non-axisymmetric term^and there- 



fore L z would be conserved ( Kozai 19621 ) . 
However, since the secondary SMBH has 
non-zero eccentricity, the perturbation has 
non-axisymmetric term, resulting in a tri- 
axial potential field. Thus, the orbits of 



many stars become chao tic {j^^ul koOO: 



Levin. Wu k. Thommesll2005l ). 

This chaotic change of the orbits of field 
stars means that they approach the secondary 
SMBH when their total angular momentum 
becomes small enough. One can assume that 
stars on both prograde and retrograde orbits 
will approach the secondary SMBH in roughly 
equal chances. The cross section of the large 
change of the binding energy, which would re- 
sult in the ejection of a star, is larger for pro- 
grade orbits, because the typical relative ve- 
locity is smaller for prograde orbits and the 
effect of gravitational focusing is larger. 

Figure 4 shows the cumulative mass of 
ejected stars in units of the mass of secondary 
SMBH. We can see that prograde stars are 
more likely to be ejected. In this calculation, 
the difference between the cumulative mass of 
ejected stars on prograde and retrograde or- 
bits was comparable to the mass of the sec- 
ondary SMBH. For T > 200Myrs, the ejec- 
tion rates of prograde and retrograde stars are 
rather similar, but that is simply because the 
eccentricity of the secondary SMBH has be- 
come large enough and therefore the relative 
difference in the cross section terms has be- 
come small. 

Figure 5 shows the evolution of the aver- 
age value of L z of stars which are bound to the 
primary SMBH at lOOMyrs and are ejected by 
200Myrs. The net change of L z of stars which 
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were initially in prograde orbit is small, while 
that of stars in retrograde orbits is large. Thus 
the primary mechanism of angular momentum 
transfer is the transfer of angular momentum 
to retrograde stars through secular perturba- 
tion and removal of them through ejection. 

To summarize, we found a quite efficient 
mechanism that causes the increase of the ec- 
centricity of unequal-mass SMBH binaries. It 
is the combination of the non-axisymmetric 
potential field of the secondary SMBH, which 
randomizes the angular momentum of field 
stars, and the selective ejection of stars with 
prograde orbits through close encounters with 
the secondary SMBH. Since either effects de- 
pends on masses of field particles, our pro- 
posed mechanism does not depend on N. 

This mechanism is effective at least for 
SMBH binaries with fairly large mass ra- 
tios (larger than 10:1). There is no other 
known mechanism that can make the eccen- 
tricity of a SMBH binary close to unity. Both 
the gravitational wave radiation and gas dy- 
namical effect reduce the eccentricity, and in- 
teraction with field stars, if we neglect the 
mechanism we described in t his paper, woul d 
not change the eccentricity ( Quinlan 19961 ). 
If the parent galaxy is globally nonspheri- 
cal, the SMBH binary could become eccen- 
tric because of the Kozai mechanism, but it 
would be hard to reach very high eccentric- 
ity values since the z-component would be 
conserved if the parent galaxy is axisymmet- 
ric. If the parent galaxy is triaxial, we could 
expect to see a significant effec t on the evo- 



lutio n of the semi-major axis (jBerczik et al 



2006h . but its effect on eccentricity is cur- 
rently unknown. Furthermore, it is difficult to 
maintain triaxi ality given the presence of the 
central SMBH dPfenniger fc de Zeeuw 1 1 19891: 



Hozumi fc Hernquistlliooi ) 



4. Discussion and Summary 

4.1. Equilibrium Value of Eccen- 
tricity 

If the mass ratio between the secondary 
SMBH and the field stars is infinite, in other 
words, in the continuous limit, we can ex- 
pect that the eccentricity of the SMBH binary 
reaches unity. In practice, since the mass and 
the number of field stars are both finite, they 
give random perturbations to the angular mo- 
mentum of SMBH binary. If we assume that 
there is thermal equilibrium between tangen- 
tial velocity of the secondary SMBH and ran- 
dom velocity of stars around the secondary, 
we can estimate the specific angular momen- 
tum Ls and the eccentricity e of the secondary 
SMBH as 



I s "~ W T77 ' Lfs> ' 



and 



1 - e 



Lt 



< L FS >) M s ' 



(2) 



(3) 



where < Lps > is the mean value of the 
specific angular momentum of the field star 
around the secondary SMBH and Ms and 
Mps are masses of the secondary and the field 
star, respectively. 

Our results in Figure 2 is consistent with 
these theoretical estimate. The eccentricity of 
the real SMBH binary can become much close 
to unity, since the mass ratio to the field star 
is much larger. 

4.2. Merging Timescale 

The timescale of the evolution of the eccen- 
tricity, T e = (1 — e)/|e|, due to our proposed 
mechanism is much smaller than the Hubble 
time. 

The timescale of the change in the angular 
momentum of the field stars is proportional to 
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M p /M s , where M p and M s are the masses of 
the primary and secondary SMBHs, respec- 
tively. The rate at which angular momen- 
tum is removed through ejection of stars is 
oc (M s /M p ) 2 , since it is proportional to the 
interaction cross section. Therefore, if M p /M s 
is large, the timescale of the change in the an- 
gular momentum is shorter than that of the 
removal of the field stars. Thus T e is deter- 
mined by the removal timescale. 

The loss rate of the angular momentum of 
the secondary SMBH is proportional to pV s a, 
where p is the density at the radius of the 
secondary, V s and a are a Kepler velocity and 
cross section of the secondary. Thus, T e is pro- 
portional to M s /(pV s cr). From our numerical 
result, we can estimate the timescale as 



s. / M v \ 3/2 / M s 
T P ~ 1.4 x 10 8 I ^ p I I - 



^ -3/2 , 

a \ ' l p 



yrs (4) 



lOpcy \W0M Q /pc 3 

o / M v \ 3/2 / M s \ 

2 0x10 (iot^J Uwfe) 

( a stall \ 3// " 
V 4pc J 



yrs, 



-2 



(5) 



where a s t a ii is a semi-major axis of a stalled 
binary and assumed a sta ii oc (Mg/p) 1 / 3 . 

Since the majority of galaxy- galaxy merger 
events are minor mergers with mass ratios 
of around 10:1 or more, the typical mass 
ratio of the SMBH binary is also around 
10:1, if we assume the linear relationship 
between the SMBH mass and the mass of 
the stellar spheroid (Magorrian et al. 19981 : 
Marconi fc Hunt! 120031 ) . Therefore, we can 
conclude that our mechanism is effective 
for most of SMBH binaries, and they can 
merge through interactions with field stars 
and gravitational wave radiation. Our result 
finally has given a de finite solution to the 
"last-parsec problem" ( Begelman et al. 1980l : 



Makino Funatoll2004l : iBerczik et all 120051 ) 



4.3. Observation of SMBH Binary 

The mechanism we found naturally ex- 
plains why the OJ287 system is eccentric. 
Most likely, this mechanism resulted in the or- 
bit of secondary with the periastron distance 
similar to the present value, but with much 
larger semi- major axis and eccentricity, and 
gravitational wave radiation have driven the 
semi-major axis to its present value. 

The state of a recently found can didate 
SMBH binary (jBoroson fc Lauerll2009l ) might 
also be explained by our mechanism. This 
system is a binary with the mass ratio of 
around 30:1. Interaction with field stars can- 
not drive its semi- major axis to the value 
suggested by the observed relative velocity 
of 3000km s -1 , if we assume a circular or- 
bit. However, our mechanism implies that all 
unequal-mass SMBH binaries with mass ra- 
tio larger than 10:1 become highly eccentric. 
Thus, they all have very large relative veloci- 
ties at the pericenter. Therefore the observed 
relative velocity is only weakly coupled to the 
semi- major axis. 

Once the binary becomes highly eccentric, 
the gravitational wave should be observed at 
each pericenter passage. The amplitude is 
given by 



h ~ 10 



-15 



( M P \ ( 



M, 



Vio lo M y Vio 8 M y uoooau 



(6) 



Vig pc ; 



The frequency / is given by 



/~io- 



M p + M s 
1O 1O M 



0.5 



1000AU 



-1.5 



Hz, 
(7) 

where r p and R are the pericenter distance of 
the secondary SMBH from the primary and 
the distance of the binary from us. Unfortu- 
nately, the frequency is too low to be observed 
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even with space missions like LISA. However, 
since the amplitude is very large, other tech- 
niques such as Doppler tracking of interplane- 
tary spacecraft or pulsar timing could be used. 

4.4. Summary 

In this letter, we present a stellar-dynamical 
mechanism that drives the increase of the ec- 
centricity of an SMBH binary. There are two 
key processes involved. The first one is the 
Kozai mechanism under non-axisymmetric 
potential, which effectively randomizes the 
angular momenta of surrounding stars. The 
other is the selective ejection of stars with 
prograde orbits. Through these two mecha- 
nisms, field stars extract the orbital angular 
momentum of the SMBH binary. Our pro- 
posed mechanism causes the increase of the 
eccentricity of most of SMBH binaries, result- 
ing in the rapid merger through gravitational 
wave radiation. 
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Fig. 2. — Evolution of eccentricity e (1 — e 
is shown) for all runs. Top and bottom 
panels show the results of simulations with 
N of 32768 and 16384, respectively. 
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Fig. 3. — Time evolution of x (top panel), 
y (middle panel) and z (bottom panel) 
components of the specific angular momen- 
tum of a star (in units of specific angu- 
lar momentum of the secondary SMBH at 
T = 0). 
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Fig. 4. — Cumulative mass of the ejected 
stars in units of the mass of the secondary 
SMBH. Solid (Dashed) curve shows the 
number for stars in prograde (retrograde) 
orbit at the time of being ejected from the 
SMBH binary. 
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Fig. 5. — Evolution of the average value 
of L z of stars, which are bound to the pri- 
mary SMBH at lOOMyrs and are ejected by 
200Myrs, in units of specific angular mo- 
mentum of the secondary SMBH at T — 
lOOMyrs. Solid and Dashed curves show 
the average value of initially prograde and 
retrograde particles, respectively. 
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